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SUMMARY/OVERVIEW: 

This  project  is  focused  on  the  study  of  flows  inside  diesel  injector  flow  passages  with 
emphasis  on  resolution  of  time-dependent  cavitation  regions  within  the  device.  During  the  past 
year,  the  turbulent  axisymmetric  model  has  undergone  substantial  validation.  In  addition,  an 
improved  inflow  boundary  condition  has  been  developed  to  increase  accuracy  in  predicting 
discharge  characteristics  of  these  devices. 

TECHNICAL  DISCUSSION 

A drilled  orifice  is  widely  used  to  provide  a simple  solution  for  atomizing  liquids  and  the 
presence  of  cavitation  inside  such  a nozzle  has  received  much  attention  in  the  past  due  to  its 
broad  applications.  Recent  experiments  have  shown  that  the  turbulence  in  the  nozzle  hole 
resulting  from  cavitation  is  a mechanism  that  promotes  atomization.  Gopalan  and  Katz1'1 
observed  that  the  unsteady  cavity  collapse  involves  substantial  increases  in  turbulence  intensity, 
and  momentum  and  displacement  thicknesses  in  the  boundary  layer.  They  also  showed  that  the 
collapse  of  bubbles  is  the  dominant  source  of  vorticity  downstream  of  a cavity.  The  present 
calculations  add  k - co  turbulence  model  to  Chen  and  Heister’ s121  homogeneous  fluid  model  to 
simulate  the  turbulent  cavitating  flow  in  a nozzle  hole.  We  chose  to  use  the  k - co  model  because 
it  has  better  performance  than  other  turbulence  models  in  adverse  pressure  gradient  flows  which 
are  common  cases  in  cavitating  flows. 


Figure  1 Schematic  representation  of  a sink  approximation  of  inflow  veloctiy. 
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In  flow  velocity  boundary  conditions  are  obtained  by  placing  an  artificial  sink  at  the 
origin  as  shown  in  Fig.  1 . For  a 2-D  problem  the  velocity  at  the  inflow  boundary  is  calculated  by 
ujn  = -(A/27tx)cos9,  Vjn  = -(A/27ir)sin9  . The  strength  of  the  sink  is  updated  during  each 
time  step  by  the  conservation  of  mass  flow  rate  through  the  nozzle  passage.  The  only  difference 
between  an  axisymmetric  problem  and  a 2-D  problem  at  this  point  is  that  a three-dimensional 
sink  is  utilized  for  the  former.  A few  axisymmetric  runs  were  made  to  compare  with  the 
discharge  coefficient,  Cd,  measurements  by  Nurick[3)  on  a circular  orifice.  Fig.  2 shows  the 


K 


Figure  2 Discharge  coefficient  Cd  comparison  with  experimental  results;  L/D  = 6 , 

D = 3.18  mm . 


Figure  3 Comparison  of  velocity  profiles  at  exit  for  laminar  and  turbulent  solutions  L/D  = 6 , 

D - 3.18  mm  . 
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results  of  a comparison  on  an  orifice  of  L/D  = 6,  D = 3.18  mm,  under  a back  pressure  of 
P2  =13.8  psi . As  shown  in  Fig.  2,  the  sink  inflow  velocity  treatment  (Lam  2)  greatly  improves 
the  prediction  of  discharge  coefficient.  The  turbulence  model  shows  a further  improvement  and 
gives  results  somewhat  closer  to  the  experimental  data.  The  differences  between  the  turbulence 
model  and  laminar  predictions  on  Cd  might  be  explained  by  Fig.  3 in  which  the  exit  velocity 
profiles  are  plotted  for  both  laminar  and  turbulent  calculations.  The  turbulent  velocity  profile  is 
fuller  than  the  laminar  one.  Due  to  this  feature,  the  turbulence  model  yields  a larger  Cd  by  1.6%. 


Figure  4 Side-View  overlay  of  cavitation  field  and  predicted  pseudo-density  contours.  The 
upper  image  provides  a comparison  with  a laminar  flow  calculation  while  the  lower  image 
depicts  the  improved  results  with  a turbulent  simulation. 

Fig.  4 shows  density  contours,  which  denotes  the  cavity  region,  obtained  from  laminar 
(the  upper  one)  and  turbulent  (the  lower  one)  calculations,  overlaying  one  photographic  snapshot 
by  Henry. Although  the  laminar  simulation  results  in  an  overall  cavitation  extent  consistent 
with  experiment,  it  indicates  two  separate  regions  of  cavitation.  The  turbulent  model  improves 
on  this  point  by  generating  a single  cavitation  region  which  appears  to  be  quite  consistent  with 
experimental  results  both  in  axial  and  cross-stream  directions. 

Fig.  5 shows  velocity  profiles  inside  the  location  of  the  cavity  region  for  cavitating  and 
non-cavitating  conditions  in  a 2-D  slot.  The  flow  inside  the  cavitation  region  in  general  is  slower 
under  cavitating  conditions  than  under  non-cavitating  conditions.  Especially  in  the  middle  of  the 
cavity(x  = 2)  a strong  reverse  flow  occurs  between  the  wall  and  cavity.  At  the  end  of  the 
cavitation  region  (x  = 3)  the  velocity  profiles  approach  the  same  shape  for  both  cavitating  and 
non-cavitating  conditions.  Cavitation  has  a significant  effect  on  the  boundary  layer  development 
downstream.  Although  the  turbulence  model  predicts  the  boundary  layer  thickness  is  almost 
identical  for  both  cavitating  and  non-cavitating  conditions,  it  produces  increased  displacement 
and  momentum  thicknesses  for  cavitating  conditions.  Other  researchers^  have  observed  similar 
behavior  experimentally. 
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Figure  5 Velocity  profiles  under  catitating  and  non-cavitating  conditions,  Re  = 23780 , K = 1.4. 


CURRENT  EFFORTS 

Most  of  the  above  results  were  published  in.[5]  Current  efforts  are  aimed  at  completing 
the  work  to  be  done  under  this  grant.  Turbulence  model  calculations  are  being  done  of  the 
nozzle  flow  experiments  of  Katz.*'1  To  date  the  comparison  shows  a discrepancy  between  the 
turbulence  model  computations  and  the  experimental  results.  Much  of  the  difficulty  seems  to  lie 
with  the  turbulence  model’s  ability  to  capture  single  phase  separated  flow.  In  addition  work  is 
progressing  on  examining  three-dimensional  effects  due  to  a cross  flow  velocity  at  the  inlet.  The 
on-going  work  should  be  completed  near  the  end  of  the  summer. 
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